Mean parameter model for the Pekar-Frohlich polaron 
in a multilayered heterostructure 

M. A. Smondyrev a )0, B. Gerlach 6 ), M. 0. Dzcro a )t 
a ) Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research 
141980 Dubna, Moscow Region, Russia 
h ) Institut fur Physik, Universitdt Dortmund, ^4.221 Dortmund, Germany 

(February 6, 2008) 



O 
O 
O 

(N 

C3 
33 
>—i 

O 

(N 

<N 
> 
On 
O 
O 

m 
o 

O 
O 



C3 

o 
o 



The polaron energy and the effective mass are calculated 
for an electron confined in a finite quantum well constructed of 
GaAs I ' Al x Ga\- x As layers. To simplify the study we suggest a 
model in which parameters of a medium are averaged over the 
ground-state wave function. The rectangular and the Rosen- 
Morse potential are used as examples. 

To describe the confined electron properties explicitly to 
the second order of perturbations in powers of the electron- 
phonon coupling constant we use the exact energy-dependent 
Green function for the Rosen-Morse confining potential. In 
the case of the rectangular potential, the sum over all inter- 
mediate virtual states is calculated. The comparison is made 
with the often used leading term approximation when only 
the ground-state is taken into account as a virtual state. It 
is shown that the results are quite different, so the incorpora- 
tion of all virtual states and especially those of the continuous 
spectrum is essential. 

Our model reproduces the correct three-dimensional 
asymptotics at both small and large widths. We obtained a 
rather monotonous behavior of the polaron energy as a func- 
tion of the confining potential width and found a peak of 
the effective mass. The comparison is made with theoreti- 
cal results by other authors. We found that our model gives 
practically the same (or very close) results as the explicit cal- 
culations for potential widths L > IOA. 

PACS numbers: 73.20.Dx + 71.38.+i 



I. INTRODUCTION 

Quasi-two-dimensional (2D) systems have attracted a 
lot of attention during the last decade because of their 
practical realization. If a heterostructure is made of polar 
materials such as GaAs, InAs etc., the electron-phonon 
interaction modifies the properties of the electron con- 
fined to a 2D-structure resulting in a shift of the binding 
energy and the effective band mass. 

The polaron effects in a 2D electron gas have exten- 
sively been studied. At earlier stages the attention was 
paid to the properties of a polaron confined to an infinite 
thin 2D layerm3. The binding energy and the effective 
mass were calculated for a GaAs/ Al~G.ai- x As infinitely 
deep quantum well of a finite lengtbcm. In these papers 
only the interaction with the bulk LO-phonon mode has 
been taken into account. Actually, LO-phonon modes are 
modified in a 2D layer (the so-called confined slab LO- 



phonon modes). Besides, there exist interface optical- 
phonon modes as wglL-as half-space LO-phonon modes 
in a barrier materialuiill For the review of these modes 
(also in complicated multi-layer-structures) see the book 
by Pokatiloy— .Epmin and BcriltJ and also more recent 
publicationstHrli^l of this group. The influence of the men- 
tioned modes on polarons were studied in Refs. [hi] |l7| 

While different phonon modes were studied in details, 
the quantum well potential was supposed to be infinitely 
deep in the cited papers. On the other hand, the proper- 
ties of the system would be quite different if a confining 
potential had a finite depth. Indeed, for an infinitely deep 
confining potential the binding energy is the monotonous 
function of the potential width which varies between lim- 
iting values E^y = ai n tkJi n for the three-dimensional 

(3D) space and E^ = (tt/2) OLiuhuJiu, where o^ n is the 
standard Frohlich electron-phonon coupling constant and 
Ldi n is the LO-phonons frequency for the quantum well 
material. If a particle is confined to a finite potential 
well, the limiting value of the binding energy should be 
the same at large width of the well. But when the width 
becomes too small, the energy level approaches the edge 
of the well, so that effectively the particle is spread over 
the 3D space. Thus, the limiting value of the binding 
energy should coincide with the 3D limiting value rather 
than with the 2D one. This means, the binding energy 



takes the value E. 



(out) 
3D 



Tiujmit at small widths where 



the parameters ui ou t and a ou t are now related to the bar- 
rier material. The binding energy evidently has a peak at 
some intermediate value of the width if -E^i) < E$p . 
If this is not the case, the existence of the peak should 
be checked in more detail. 

Different rectangular quantum wells of a finite heig, 
have been investigated by Hai, Peeters and Devreesaifrl 
and Shi, Zhu et al.E] in the scope of the second order per- 
turbation theory in powers of the electron-phonon cou- 
pling constant a with all phonon modes being incorpo- 
rated. Peaks of the phonon induced energy shift and the 
polaron effective mass were found for some values of the 
confining potential widths. 

In principle, the same approach can be used while deal- 
ing with a quantum well constructed of layers of differ- 
ent materials. But the problem becomes then too com- 
plicated because one has to take into account interface 
phonon modes at each frontier of different materials as 
well as quantized phonon modes inside each of the layers. 



1 



The main goal of the present paper is to formulate a sim- 
plified model to take these effects into account and to deal 
with the effective confining potential and only one bulk 
phonon mode. We calculate polaron characteristics for 
the same rectangular quantum well as in Refs. p"8|-p0| to 
compare the results. Another example is given of a quan- 
tum well of the finite depth for which the second-order 
correction due to the electron-phonon interaction can be 
calculated explicitly. Namely, we take the Rosen-Morse 
potential to confine electrons to a 2D-multilayered het- 
erostructure and calculate the shift of the ground-state 
energy and the effective mass perturbatively, that is, in 
the weak-coupling limit. In contrast with the rectangular 
potential we should not worry about the correct including 
of all virtual states because the Green function is known 
analytically for this system. 



II. FORMULATION OF THE MODEL 

Let us consider a quantum well in the z direction con- 
structed of the xy plane layers of GaAs / Al x Ga\^ x As. 
That is, the AlAs mole fraction x depends on the co- 
ordinate z: x — x{z). The energy gap between different 
materials forms the confining potential V(z) which serves 
us as the main entity. Given the potential V(z), one can 
find the corresponding mole fractions x(z) and a depen- 
dence on z of any of the medium parameters (such as 
the electron band mass m(z), phonon frequencies lu(z), 
dielectric constants £o(z), £00(2;), Frohlich coupling con- 
stants a(z), etc.). 

To avoid difficulties with mass mismatch in different 
layers we suggest to use a mean band mass to which is 
common for all layers. Then we start with the following 
expression for the electronic part of the Hamiltonian: 



H e l = #ei,|| + Hel,±, 



H, 



el.\\ 



2m 



H, 



(2.1) 



where the electron mean band mass is defined by the 
relation 



1= ( dz \M^L 

to J m(z) 



(2.2) 



and the ground state wave function ipi (z) for the electron 
motion in z direction is a solution to the Schrodinger 
equation 



(2.3) 



with Ei being a ground state energy. As the wave func- 
tion also depends on the mean band mass to, the l atte r 
can be found as a self-consistent solution of Eqs. (2.2), 

s 

In a similar way we define the free LO-phonon Hamil- 
tonian 



(2.4) 



where at (ag) are creation (annihilation) operators of a 

phonon with a wave vector k, and mean frequency w LO 
can be found from the expression 



w LO = / dzuj(z) \ipi(z)f 



(2.5) 



Evidently, we have to address why the free phonon 
Hamiltonian is averaged with respect to the electron wave 
function. Our motivation is based on the fact that we 
are going to apply our model to calculate polaron ef- 
fects. That is, our effective phonons will be considered 
only in a cloud around the electron, and the properties 
of this cloud depend on the electron position. So, in our 
model the effective phonons replace numerous phonon 
modes whose frequencies depend on the coordinate z of 
the electron. 

Finally, we accept the conventional form of the Hamil- 
tonian describing the interaction of the electron with ef- 
fective phonons: 



Hint. — 



k 



a l y* e -ik-f 
k k 



(2.6) 



where the Fourier transforms of the electron-phonon in- 
teraction potential are specified as follows: 



Ana 



Vk 2 V 2muj L c 



1/2 



(2.7) 



Here the use is made of a mean Frohlich coupling constant 
a which can be found from the relation 



'a 



dz\Mz)\ 2 — 
W LO 



a(z) 



muj LO 
m{z)uj{z) 



1/2 



(2. 



Note t hat we define the mean parameters in Eqs. (p. 2]), 
( |2.5| ), (2.8) according to the way they enter the Hamilto- 
nian. 

Thus, we describe a complicated multilayered het- 
erostructure by the Hamiltonian 



H = H e l 



H 



ph 



Hi, 



(2.9) 



with the bulk phonon mode only which inhabits an ef- 
fective medium with mean characteristics defined above. 
The details of the heterostructure are taken into account 
by the confining potential V(z). 

Performing a unitary transformation H — > H' = 
U~ l HU with the operator 



U = exp 



5^11 4 a fe 



(2.10) 



we arrive at the Hamiltonian 
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H' = H' elAl + H eli± + H ph + H' int , (2.11a) 

2 



^,11 - 2i „ 



at V; e~ lk * 

k k 



(2.11c) 



The quantity p\\ is a c-number corresponding to the con- 
served momentum in the xy pla ne a n d th e Hamiltonians 
H e i,±,H p h are defined by Eqs. ( |2.l| ), (2.4), respectively. 

Keeping in mind the smallness of the electron-phonon 
coupling constant a for most of the materials, we cal- 
culate the second-order correction to the unperturbed 
Hamiltonian H' Q = H' el ^ + H e i t ± + H p h (note that the 
quantum-mechanical first-order correction is equal to 
zero). The unperturbed energy levels are given by the 
expression 



1 



E(p h n^,N) = —\p 



n Y,^ n k 



k 



■ E 



N, 



(2.12) 



where is the number of phonons with the wave vector 

k. The energy Em is the iV-th energy level of the one- 
dimensional system H e i y ± of Eq. ([j~l|). Here N is the cor- 
responding quantum number not necessarily a discrete 
one: it stands for both the quantum number n which 
varies from 1 to n max and the wave vector q of the con- 
tinuous spectrum states. 

The wave functions of the unperturbed Hamiltonian 
Hq are given by the direct product 



\pfn s ,N) = \n n )®\N) 



(2.13) 



of the corresponding wave functions of different terms in 
H>. 

Because of the structure of the interaction term H' int 
only intermediate states with one phonon contribute to 
the second order correction to the ground-state energy. 
The latter is then given by the expression 



A 2 £(P||) 



\G(N, fc z )| 5 



N,k 



En + ftw L 



(p\\ -ftfcll) 2 -p*ll 



(2.14) 



2m 



E x 



where 



G(N, k z 



dz ij)N{z)il)i(z)e 



ik z z 



(2.15) 



and tPn{z) are the wave functions of the Hamiltonian 
H e i : ± in Eq. (2.1). The concrete application of these 
formulae is given in the following section. 



III. RECTANGULAR POTENTIAL 

A. Medium mean characteristics 

As an example we now consider the rectangular con- 
fining potential 



and 



V(z) 



m{z) 



0, \z\<L/2 
Vo, \z\>L/2 



m in , \z\ < L/2 
m out , \z\ > L/2 



(3.1) 



(3.2) 



with rrii n (m ou t) being the electron band masses in the 
well (barrier) material, respectively. For concreteness 
we assume GaAs to be the quantum well material and 
Al x Ga\- x As to be the barrier material. 

Symmetrical wave functions of the discrete spectrum in 
the rectangular quantum well with the mean band mass 
m take the form 



where 



cos q n z, 

cos^e-PnCI*!-^), 



Pn = 



\z\ < L/2 
\z\>L/2, 

2mV 



and the normalization constant is given by 



2p n 



PnL 



(3.3) 



(3.4) 



(3.5) 



Antisymmetrical wave functions of the discrete spec- 
trum take the form 



ipa,; 



I sm q n z, 

N s ,n \ / \ ■ QnL -p„(\z\-L/2) 

sgn(z) sm — — e Pn ^^' u i*)^ 



\z\ < L/2 
\z\ > L/2 
(3.6) 



with the same normalization constant given by Eq. ( |3. 5] ). 

The total number n max of the discrete energy levels is 
given by the expression 



1 



(3.7) 



where [A] is an integer part of a number A. The expres- 
sion for the discrete energy levels reads as follows: 



— atan < 



n= 1,2,.... (3.8) 



Energies with odd (even) n correspond to the symmetri- 
cal (antisymmetrical) wave functions. 
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The energy E q = h 2 q 2 /2m of the continuous spectrum 
state depends on the wave vector q. The corresponding 
symmetrical wave functions are as follows: 



N, 



s.q 



s.q 



p cos qz, 



qL 



\z\<L/2, 



pcos — cosp(|z| — L/2) — 
gsin^-sinp(| 2 |-L/2), \z\ > L/2, 



where 



P : 



(3.9) 



(3.10) 



and L z is the (infinite) size of the system in the z direc- 
tion. The normalization constant is given by the expres- 
sion 



2V, 



\ 2 2 , 2 • 2 ?^ 

^ cos^ — + 9 sin — 



(3.11) 



The antisymmetrical wave functions are as follows: 
' Psinq\z\, \z\<L/2, 



1pa,q = 



N a ^ q sgn(z) 



psin %- cosp(|z| — L/2)+ 
q cos ^- sin p(\z\ - L/2), 



\z\ > L/2, 



(3.12) 



To finish this subsection, we note that the exact energy 
levels in the rectangular potential with different masses 
rriin and m ou t calculated for the GaAs / Al x Ga\^ x As het- 
erostructure practically coincide with the levels obtained 
with the electron mean band mass m. To obtain an in- 
ner criterion of the validity of the anzatz concerning the 
mean band mass we notice that the particle being on low- 
est energy levels is located mostly inside the well which 
means that its band mass is almost coincide with nii n . 
One can await the largest discrepancy for a level near 
the potential edge. The n-th discrete level appears at 
the width L = L n av \ where 



n(n — 1) 



(3.16) 



and the analogous width for the exact solution reads as 
follows: 



4«> = An 1) 



h 



y/2m in Vo 



Thus, the ratio 

r (av) 



j- {ex) 



(1- Wi, 



m out 



(3.17) 



(3.18) 



can serve us as the numerical criterion of the validity of 
the anzatz. The largest discrepancy happens at n = 2 
and in this case Eqs. (|3.4|), (3.8) (3.15) lead to the 
following expression: 



= J 0.844 + 0.156- 



(3.19) 



where the normalization constant is given by the expres- 
sion 



N a ,n = 



\ 9 • 2 QL i 2 2 9-^ 

^ p z sin — + q cos — 



(3.13) 



The electron mean band mass is defined as 

1 _ W m | W out 
m m m rtiout 



W in m out + (1 - W in )mi 



(3.14) 



where Wi n and W ou t — 1 — Wi n are probabilities to find 
the electron inside (outside) the quantum well. The ex- 
pression for Wi n follows from Eq. (|3.3[) 



(3.15) 



L/2 

Wi n = 2^,1 J dz cos 2 q x z = 1 - -| 



{qi/q-maxf 



PiL/2 



where qi is a solution to Eq. ( p.8[ ) for the ground-state 
(n = l). 



Note that numerical coefficients here do not depend on 
the material parameters. For the GaAs / AIq^GclqiAs 
quantum well we have mi n /m out w 0.7 and the dis- 
crepancy is about 2%; in the worst possible case when 
Triin/rriout <C 1 the discrepancy is still not large: 
100%v'O44« 8%. 



B. Electron-phonon correction to the polaron energy 
and the effective mass 



Summation over the wave vector k in Eq. (2.14) can 
be reduced to integration in a conventional way 



Ei^i 2 (■■•) 



V 



(27T) 



, , dk\Vj:\ 2 (...) 



dku dk z 
2muj LO 2tt 2 J B + kl 



(...). 



(3.20) 



Then, the integration over fc|| in Eq. ( |2.14j ) can be per- 
formed explicitly. As we are interested in corrections to 
the ground-state energy and the effective mass J7i e // f» 
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m + A 2 m of the polaron motion in the xy plane, we 
expand A2-E(pj|) as A 2 E — ^ 2 2 p| in powers of the con- 
served momentum pj| . Doing this the use is made of the 
integral 



fc 2 -& 2 



(fc 2 + fc 2 )[fc | 2 -2fc r p||/?i + 6 2 ] 

P|\ 2 fc| - b 4 - 2fcgb 2 ln(fc 2 /fr 2 ) 
?i J ^ b 2 (k 2 -b 2 ) 3 



(3.21) 



As the next step we use dimension less "pola ronic" 
units performing the scaling k z —> k z y / 2muJi j0 /h, z —> 
zy/h/2moj LO and using the notation 



2mu! LO E N 
l=L\l - , £N = t • 



(3.22) 



In these units the correction to the ground-state energy 
takes the form 



A 2 E 



00 

n i 



k 2 - h 2 
K z U N 



{\G s (N,k z )\ 2 + \G a (N, k z )\ 2 ), 

(3.23) 



where 



b N = Vew + l-ei- (3.24) 

The correction to the effective mass reads as follows: 
A2TO 



DC 



b 4 Ar -2fc 2 b 2 r ln(fc 2 /6 2 v 



(|G s (7V,fc z )| 2 + |G a (iV,fc z )| 2 



(3.25) 



Quantities G 3 (N,k z ) in Eqs. ( ^23| ), C25h are given in 
dimensionless units by the same Eq. (2.15); the indices 
(a)s are related to (a nti)symmctrical wave functions be- 
ing used in Eq. (2.15): 



G s (N,k z ) — 2 J dztp Si N(z)tp s ,i{z) cosfc z z, 



00 

G a {N,k z ) = 2 J dz<ip a , N {z)iJ Sll (z) sinfc z z, (3.26) 


Evidently, the replacement L — > I should be done in the 
definition of the wave functions and their normalization 



constants; in addition L z — > l z in Eqs. (p-9), (3.12) as 



well as in Eq. (|3.8|) for the energy levels of the discrete 
spectrum. Eq. (|3.4j) now reads as follows: 



Pn = \/vo- g 2 , v 



Vo 



Vq. 



Eq. (3/7) takes the form 

Tlm.nm 1 H 



v l 



(3.27) 



(3.28) 



The relation of dimensionless energies of the discrete and 
continuous spectra with subsequent wave vectors takes 



the form e n 



2 



q . All the changes mentioned 



should also be done in Eq. ( |3.15| ). 

The final note of this section concerns summation over 



N in Eqs. ( |3.23|) , ( p.25| ) 



n=l 



■) + lim — dp(. 

l z —*-C>Q Z7T 



(3.29) 



The replacing of the sum over the continuous spectrum 
by the integration over the wave vector p follows from 
Eqs. ( 3~9| , 3.12j ) in the limit L z — > oo. The wave vectors 
q and p are related to each other because of Eq. ( 3.10| ) 
which now takes the form q = p 2 + vo . Note also that 
only G S (N, k z ) (G a (N, k z )) has to be taken into account 
for odd (even) n in the sum over the discrete quantum 
number n. 

The numerical results obtained are plotted in Fig. |l| 
for A2-E and in Fig. ||a for Aimj-m. Because the mean 
effective mass m depends on the potential width we also 
plotted in Fig. ||b the ratio of the mass shift to those in 
the well material, that is, the ratio 

x A2m a am 11 qr^ 

o 2 m = , A 2 m in = m m ——. (3.30) 

A 2 m m 6 

The discussion of the numerical results is given in the 
last section. 



IV. ROSEN-MORSE POTENTIAL 

A. Energy-dependent Green function 

In this section we present another example — a mul- 
tilayered heterostructure described by a confining poten- 
tial V(z) which is chosen in the form of the Rosen-Morse 
potential 



V(z) = Vbtanh 2 



Vn = 



J RM 



2mL 2 RM 



k(k + 1). 



(4.1) 



where Lrm is the parameter close to the half-width of 
the Rosen- Morse quantum well and k is the dimensionless 
parameter to govern the strength of the potential. 
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The summation (2.14) over the quantum number N 
can be represented through the Green function which 
is known analytically for the Rosen-Morse potential. 
Namely, the second-order correction to the ground-state 
energy can be written in the form 



A 2 E = -%uj. 



l>RM 

V2 



dkn 



dz a / dz b e - fc ii !*•-*»! 



ip$(z a )i>i(zb) G(z a ,z b ;E), 



(4.2) 



where we made a scaling z — > zLrm-, k — > k/LuM to use 
dimensionless variables z, k and integrated over k z and 
angles of feii . The dimensionless parameter 



Irm = Lrm 



(4.3) 



is the width of the confining potential in polaronic units 
while the potential strength can now be written as fol- 
lows: 



V = hu L 



k(k + 1) 

~ 2I 2 ' 

M RM 



(4.4) 



The quantity G{z ai z bl E) is t he Green function of the 
dimensionless Hamiltonian (|2.l[ ) which takes the form 



tti, 1 d2 k(k + 1V , 2 



(4.5) 



that is G(z a ,z b ;E) = {z a \{H'^ ± - E)~ 1 \z b ), while ^{z) 



is the ground-state wave function of the potential (4.5) 
T(k + 1/2) 1 1/2 



Mz) 



1 

cosh'" 2 



(4.6) 



The ground-state energy of the Hamiltonian fl4.q ) is given 
by 



E, = -. 



The energy E in Eq. (|4.2|) reads as follows 



k 2 

^ — l RM 



K 

2' 



(4.7) 



(4.8) 



The energy-dependent Gre cp . fu nction of the system 



can be represented in the forml 
G(z a ,z b ;E) = 

T(v + k + l)T(y - k) 



a 



i 



r 2 (^ + i) 

2 Fi [v — k, v + k + 1 ; v + 1 



(4 cosh z a cosh 25)^ 
1 — tanh z > 



1F1 v — k, v + k + 1; v + 1; 



1 + tanh z< 



(4.9) 



where z> (z<) denotes the maximum (minimum) of z Q 
and Zf,. The parameter v is defined by the relation 



v = \ -2\ E 



k(k + 1) 



= Jkf 



n 2 + 2l RM . (4.10) 



The polaron effective mass can be represented in a sim- 
ilar way 



I 

= a 



_ j dk\\ fcff / dz a / dz b e~ k * - ; 1 




^(^i(^) -^G(z a ,z b ;E). (4.11) 

To simplify numerical calculations we may replace the 
derivative with respect to E by the derivative with re- 
spect to v 



d 2 
dE 2 



v 2 dv 2 



1__9_ 

v z dv 



(4.12) 



and perform once the integration by parts. As the result, 
we arrive at the following representation equivalent to 
Eq. pi: 



Ayin 



dkn 



dz n 



dz b (1 - k\\\z a - Zb\ 



-k«\z a -z b \j* 



ipl(z a )ipi(z b ) 



ld_ 

v dv 



G(z a ,z b ;E). (4.13) 



Note that to, a, cj lo in all these formulae stand for the 
mean characteristics of the medium. The wave function 
in their definitions is given by Eq. ( |4.6| ). The numerical 
results are plotted in Fig. || and discussed in the last 
Section. 



B. Effective width 

If we decide to compare the results for the rectangu- 
lar and Rosen-Morse potentials, we have to define a pa- 
rameter which plays the role of the effective width of 
the Rosen-Morse potential. That is, this parameter (for 
which we use a notation L) should be close to 2Lrm 



of Eq. (4.1) being also related to the rectangular poten- 
tial. We accept the following definition: let us call the 
effective width of the Rosen-Morse potential the width 
L of the rectangular well of the same height with the 
same ground-state energy in the absence of the electron- 
phonon interaction (that is, at a = 0). The advantage 
of this definition is that while calculating the polaron 
binding energy for the Rosen-Morse and rectangular po- 
tentials we subtract the same quantity in both the cases 
and can compare only energy shifts due to the electron- 
phonon interaction. 



G 



The ground-state energy of a rectangular potential 
with the height Vq and width L is given by the relations 



2l2 



Erc = 
kL 



2m ' 



tan • 



To 



1, 



2 V E R c 
while the RM ground-state energy looks like 



E 



RM 



m ^RM 



(4.14) 



(4.15) 



and the height Vq of the potential is given by Eq. (4.1). 
With the equality Erm = Erc we arrive at the relation 
between the parameter Lrm of the Rosen-Morse poten- 
tial and its effective width defined as has been discussed: 



— = 2\/A arctgV A — 1, 



A = K 



1 

1= 2 



1 + Jl + (2L RM /L 



Here we introduce the distance scale 



L 



2mV 



(4.16) 



(4.17) 



The rel atio n to the other dimensionless parameter Irm 
of Eq. (h3) is given by 



^RM 

L 



2V Q 



kUJr 



(4.18) 



At sm all Lrm <C Lq wc obtain L ~ 2Lrm from Eq. 
( 4.16 ), that is indeed the parameter Lrm plays a role 



of the half-width of the Rosen-Morse potent ial in this 
case. When Lrm ^> La, it follows from Eq. ( 4.16] ) that 
L ~ h\JLrmLq. 

The effective width L defined in this subsection allows 
us to apply the results for the rectangular potential to 
the Rosen-Morse quantum well. The example is given 
in Fig. U where we plotted also the energy and the mass 
shifts for the rectangular potential vs. the parameter 
Lrm related to L as is described. 



V. NUMERICAL RESULTS AND DISCUSSION 

To proceed to the numerical calculations we need now 
the dependence of medium parameters on the AlAs mole 
fraction x. At first we_present the parametrization from 
the review by AdachiO: 



a(z) = 0.068 + 0.058a;, 



m(z) = m e ■ (0.0665 + 0.0835a;), 



(5.1a) 



(5.1b) 



%w(z) = (36.25 + 1.83a: + 17.12a; 2 - 5.11a; 3 ) meV, (5.1c) 



which was us 



numerical calculations by Hai, Peeters 
Here m e is the electron mass in vac- 
its band mass in the subsequent layer; 



and Devreesefc 
uum and m(z) 
the values of the electron-phonon coupling constant ct(z) 
and the LO-phonon frequency lu(z) are also related to 
this layer. 

Some comments are to the point. The expression for 
the electron band mass is nothing else but the linear in- 
terpolation between the values m = 0.0665m e for GaAs 
and m = 0.150r7i e for AlAs. As to the LO-phonon 
frequency there are two phonon modes with different 
frequencies lo^ g \z) and uj^(z) for the GaAs-like and 
AlAs-like modes in Al x Ga\- x As crystal. Experimental 
results of Ref. |24| are interpolated by the following for- 
mulae: 



huj {G \z) = (36.25 - 6.55a; + 1.79a; 2 ) mcV, 



(5.2a) 



huj iA) (z) = (44.63 + 8.78x - 3.32a; 2 ) meV. (5.2b) 

Because the exact theory of the two-phonon interaction 
in alloys where there are two-mode phonons present has 
not been reported, Adachi suggested to use the effective 

( A \ that is the 



(G) 



XOJ 



phonon frequency w = (1 — x)u> 
linear interpolation betw een these two modes. Ins ertin g 
here the expressions (|5.2|) one arriv es at the result ( 5.1c ) . 

As to the interpolation formula (5.1a) for the Frohlich 
coupling constant a, the situation seems to be a bit in- 
consistent. Indeed, a depends on the values of the static 
£o and the high-frequency dielectric constants: 



J_i_ 

hu) ^2 
116.643 



muj 



1 



1 

E0 



1 



ni 1 1 meV 



(5.3) 



Earlier measurements of Eq of GaAs have yielded widely 
different values ranging from 9.8 to 13.3 (see Ref-. |25| 
and references therein). For instance, Kartheusero re- 
ports the values Soo — 10.9 and Eq = 12.83 and hui = 
36.75 meV for GaAs. This leads to the result a — 0.068, 
which is widely known and used by many people. 

On the|Other hand, Adachi used the more recent results 
for G-svlsEZl: e = 13.18 ± 0.40 and = 10.89, and for 
AlAsM: £ = 10.06 ± 0.04 and £oo = 8JJ5 ± 0.02. This 
gives birth to his interpolation formulaeEi 



e Q = 13.18- 3.12a;, 



10.89- 2. 73x, 



(5.4a) 



(5.4b) 



Inserting formulae ( |5.1b ), (5.1c) and (5.4) into Eq. ( |5.3| ) 
Adachi declared the result a = 0.126 for AlAs. To- 
gether with the value a = 0.068 rep orted in Ref. ^6] this 
leads to the interpolation formulae ( [5. la ). The problem 
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is that both these values for a do not follow from the 
parametrizations mentioned above. 

Taking the same values for AlAs as Adachi did take 
(m = 0.150m e , Tilo — 50.09 meV, Eq — 8.16, = 
10.06) we arrive at the result a = 0.1477. Moreover, if 



one takes the same interpolation formulae (|5.4| ) at x = 
one obtains the value a — 0.0797 for GaAs. That is, 
Adachi had to obtain the formula 



a(z) = 0.0797 + 0.0680a; 



(5.5) 



as a linear interpolation between the values of a in GaAs 
and AlAs. Note, that this formulae can be presented 
in the form a(z) = 1.172(0.068 + 0.058a;). The expres- 
sion between the brackets coincide (probably occasion- 
ally) with the Adachi interpolation for mula e for a (cf. 
Eq. ( |5.1a| )). That is, the discrepancy of (5.1a) and of our 
interpolation ( |5.5| ) is about 17% and do not depend on x. 
To b e consistent we have to accept the parametrization 
(5J3) in what follows. 

For the confining potential we take the expression de- 
rived from the band-gap energy fit of Ref. ^9] and used in 
Rcf. 



V{z) = 600 • (1.155a; + 0.37a; 2 ) meV. 



(5.6) 



Thus , we use the parametrization (5.1b), (5.1c), ( |5. 4a ) 
( 5.4b| ), ( [5~5| ) and the potential ( |5.6| ) in our numerical cal- 
culations. 
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FIG. 1. Total energy shift A2E is shown vs. the rectangu- 
lar potential width L for x — 0.3. The contribution of the 
discrete (Adi S E) and continuous (A con E) spectra are pre- 
sented separately as well as the result of the leading term 
approximation (AuE). 

The results of our study for the rectangular potential 
(which is formed by a layer of GaAs / Al x Gai_ x As) are 
shown in Fig. |l| for the polaronic energy shift and in Fig. |^ 
for the polaron effective mass at the AlAs mole fraction 
x = 0.3. The contribution of the discrete and continu- 
ous spectra are plotted separately for this potential. In 
Fig. Ha the relative mass shift ^m/m is shown where 



the mean mass m also depends on the potential width L. 
Thus, the ratio S^m = / ' Aimi n of the mass shifts in 
the potential and in GaAs is presented also in Fig. ||b for 
the same AlAs mole fraction. Evidently, the asymptotics 
of this curve is equal to the unity at large L and to the 
ratio m out a ou t/m m a in at L — ► 0. 



GaAs/AI 03 Ga 07 As 
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FIG. 2. The relative shifts A 2 m/m 

and 82m — Aimj Airriin of the effective polaron mass for the 
rectangular potential at x = 0.3. Contributions of the dis- 
crete and continuous spectra are shown as well as the result 
of the leading term approximation. 

We may conclude that the continuous spectrum domi- 
nates at small potential widths. At large widths its con- 
tribution could also be significant although it is smaller 
than the contribution of the discrete spectrum (especially 
in deep potential wells). We also confirm the conclusion 
of the preceding papers that the leading term approxi- 
mation is not adequate to describe this system and leads 
to wrong asymptotics at both small and large potential 
widths (see the dashed lines in Figs [|, [2]) . 

An example of a multilayered heterostructure is pre- 
sented. The results for the energy and the effective mass 
for the polaron in the Rosen-Morse potential well are 
shown in Fig. ||. For the numerical calculations we fix 
the value Vo = 227.9 meV in Eqs. flO]) , (pLl]) which cor- 
responds to the limiting mole fraction at large distances 
Xrxj = lin^^oo x(z) — 0.3. Thus, we obtain the depen- 
dence of the mole fraction x on the coordinate z: 
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600 • (1.155a: + 0.37x 2 ) = 227.9 tanh 2 z. (5.7) 

Now Eqs. ( |5. lb ), (5.1c), allow one to define the 

dependence of parameters on the coordinate z and to 
calculate the mean characteristics of the heterostructure. 

The calculations were completely different in compar- 
ison with the rectangular potential: instead of the direct 
summation over all intermediate states we used the ana- 
lytical expression for the Green function of the Rosen- 
Morse potential. The results obtained demonstrate a 
similar behavior which is also close numerically to the 
results for the rectangular potential. The polaronic en- 
ergy and mass shifts for the rectangular quantum well 
are also plotted here (dashed line) vs. the Rosen-Morse 
width Lrm obtained from L as is described above. We 
see that both the energies almost coincide, which gives 
the opportunity to approximate different quantum wells 
by the rectangular potential. The discrepancy in the ef- 
fective mass is larger but not so crucial. This serves also 
as an additional internal criterion of the validity of our 
calculations. 
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FIG. 3. Polaron energy and effective mass shifts for the 
Rosen-Morse potential (solid curves). The dashed lines 
present results for the rectangular potential at x = 0.3 as 
functions of 2Lrm recalculated from the width L as is de- 
scribed in the text. 

Thus, we obtained a monotonous behavior of A2-E 
between the correct 3D limiting values cti n huii n and 
a out hio out both for the rectangular and the Rosen-Morse 
potentials (see Fig. [l] and Fig. |]a). Actually the peaks 
are "hidden" and they reveal themselves if we plot the 



dimensionless energy shift A2E / (Tilj^qo) which has the 
same 3D limit (the unity) at both small and large poten- 
tial widths. But in the "real" units (meV) the peaks are 
smoothed. 
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-Our result 

■ Hai, Peeters & Devreese 
PRB 48, 4666 (1993) 

- Hai, Peeters & Devreese 
(private communication) 
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FIG. 4. Comparison of the results of the present paper 
(solid lines) and those of Ref. |l8| (dotted line) and of Ref. |l9| 
(dashed lines) for the rectangular potential generated by the 
GaAsfAlo.3Gao.7As heterostructure. For this plot we used in 
our calculations the same parametrization (^ij) and (5.6) as 
these authors did. 

To compare our results with the calculations per- 
forme d fbr . the one-layer heterostructure we refer to the 
papersMLa where the authors took into account the con- 
tributions of different phonon modes as well as mass and 
dielectric constant mismatches in the materials of the 
barrier and the well. Note that the analytical formulas of 
Ref. [l^ contain a mistake — the wrong expression for the 
density of states. Namely, in some parts of the continuous 
spectrum contribution the integration is performed not 
over the wave vector p but over the wave vector q (that is 
dE z l\[W z (...) in the notations of that paper instead 

of the correct integration dE z j\JE z — Vq (. . .)). It is 
clear that this mistake results in lowering of the resulting 
curve for the energy, and the discrepancy is larger when 
the energy is closer to the potential edge, that is at small 
widths. This is just what we see in Fig. ^a comparing the 
result of Ref. [ll] (the curve AhpdE) with the new cal- 
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culations of the same authors (the curve AhpdE) which 
came to our knowledge when the present paper was al- 
ready submitted for the publication. 

Thus, our model does not reproduce the more compli- 
cated structure with the peak and the dip which was ob- 
tained in Ref. Some hints on the existence of peaks 
can also be seen in our plots but the maximal values 
are so close to the asymptotics that the peaks are al- 
most invisible. Probably, the dip appears because of the 
presence of several phonon modes (bulk, interface, etc.). 
At widths L > 50 A our results for the energy practi- 
cally coincide with those of Ref. |l9|. The discrepancy at 
smaller widths seems to be more crucial. But the differ- 
ence between the values in the peak and in the dip for 
the curve A H E in Fig. |a is about 0.1 eV (3%). This 
phenomena hardly can be seen experimentally and this 
discrepancy is in the limits of the accuracy of our model 
estimated above. This gives indeed a strong support to 
our model and we may conclude that the latter provides 
us with the rather accurate approximation and can be 
used for more complicated calculations in multilayered 
heterostructures. 

As to the shift of the electron band mass we found 
clear peaks for both the rectangular and the Rosen-Morse 
potentials (see Figs. ||,||). As is seen in Fig. [2] the effective 
mass shift for the polaron in the rectangular quantum 
well has a peak at L « 20 A (x — 0.3). Calculations 
show also that the larger is x the smaller is the potential 
width corresponding to the peak. For the Rosen-Morse 
potential at x max = 0.3 the peaks in the effective mass 
occur at 2Lrm ~ 20A. Note, that again the authors of 
Ref. [l^ obtained curves with peaks and dips in contrast 
with our results (see Fig. ^b). The maximal discrepancy 
for the mass is about 11% at L « 3 A which is beyond 
the region available for experiments. Our results are very 
close to those of Ref. [l^ at L > 10A and practically 
coincde with them at L > 20A. 

To compare our results with those of Ref. ^0] we 
need now another parametrization used by these authors 
(although they refer also to the paper by Adachic3). 
Namely, they took a slightly different expression for the 
confining potential: 



V(z) = 600 • (1.266a: + 0.26x 2 ) meV, 



(5.8) 



which follows from the band gap of Ref. Further- 
more, instead of the effect ive LO-phonon frequency they 
used the expression ( 5.2a ) for the energy of the Ga As- 
like phonons. The Frohlich coupling constant a wa s cal- 
culated then using also the parametrization (5.1b) and 
(5.4). Note, that these numerical calculation, as we 
found, can be approximated by the interpolation formula 



a{z) = 0.0797 + 0.0772.T + 0.0295a; 2 . 



(5.9) 



The results of the comparison are shown in Fig. |5| 
(we used in our calculations for this plot the same 
parametrization as was used in Ref. pQ) . 



The curve AcuE in Fig. ||a for the energy shift taken 
from Ref. ^ has also a small dip (qualitatively similar 
to this of Ref. [l9|). But the discrepancy between energy 
shifts is much more drastic in this case, and we have no 
explanation for this. 
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- Our result 

J.J.Shi et at, PRB 55, 4670 (1997) 
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FIG. 5. Comparison of the results of the present paper and 
those of Ref. |2Tj for the rectangular potential generated by 
the GaAs / Alo.25Gao.75As heterostructure. For this plot we 
used the parametrization of these authors as is described in 
the text. 

It is clear that at large potential width only a bulk 
phonon mode inside the quantum well contributes so 
these curves should have the same limiting value a.i n 0Ji n . 
Numerically we found a, n = 0.0797 and uii n = 36.25 
meV, so atnOJin = 2.89 meV. Moreover, the behavior of 
the curves at large L should be qualitatively and quanti- 
tatively the same which was the case when we compared 
our model with Refs. [jj|,[l9[ In contrast with our model 
and the cited results by Hai, Peeters and Devreese the 
curve AchE in Fig. [5)3 approaches the asymptotics from 
below and the subsequent mechanism remains unclear. 
On the other hand there are some reasons why the curve 
have to approach its asymptotics from above. Indeed, at 
large potential width the particle does not feel yet the 
finite height of the potential, and the energy shift takes 
the same value as in the infinitely high potential which 
is a bit larger than the free polaron energy. 

As to the opposite limit of the small width of the con- 
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fining potential, it is surprising that the asymptotic value 
is not reached even at L ~ 0.3 A, as is found in Ref. |2(i| . 
Numerically we obtained a ou t — 0.1014 in this scheme 
and oJout = 34.72 meV, so a ou tU> ou t = 3.52 meV. 

Both asymptotic values coincide with what was ob- 
tained by the authors of Ref. 20 , Looking at the behav- 
ior of the mass shift, we see that both curves coincide at 
large widths as it should be. At widths smaller than 100 
A the discrepancy becomes evident. But we may con- 
clude that something is wrong with the numerical job of 
Ref. 20 because their curve approaches the wrong limit 
at L — ► 0. Indeed, in this limit the asymptotical value 
of the plotted ratio should be equal to a ou tm out / ai n rni n . 
As it follows from our analysis of the energy shift, we 
obtained the same values for the Frohlich coupling con- 
stants^ The values for the band masses follow from 

Eq. ( 5.1b| ): rrii n = 0.0665m e and m ou t — 0.0874m e at 
x = 0.25. Then, the asymptotical value of the plotted 
ratio should be equal to 1.67, instead of 1.83 what was 
found in Ref. |2(] That is, the discrepancy is about 10% 
in this limit and we cannot explain its origin as well. 

It would be highly desirable to include a comparison 
of our results and the results by other authors with cor- 
responding experiments. To the best of our knowledge, 
no such experiments do exist at the moment. 



VI. CONCLUSIONS 

To conclude, we suggested an approximate model to 
describe a multilayered GaAs / Al x Ga\- x As heterostruc- 
ture as an effective medium with one (bulk) phonon 
mode. The fundamental entity is the confining potential 
generated by these layers which we take into account ex- 
plicitly. Then we calculate the mean characteristics of the 
electron in the effective medium (such as its band mass, 
phonon frequencies etc.) which depend on the form of 
the confining potential. With these parameters we calcu- 
lated the energy and the effective mass of a polaron con- 
fined to a quasi-2D quantum well GaAs / Al x Ga\- x As for 
different AlAs mole fractions. The calculations include 
the full energy spectrum as intermediate states. Peaks 
are found for the effective mass at some potential widths 
while the energy demonstrates rather monotonous behav- 
ior between the correct 3D-limits. Finally, some discrep- 
ancies in the interpolation formulae for the experimental 
results are discussed as well as discrepancies with the pre- 
viously obtained theoretical results. We demonstrated 
that our model gives practically the same (or very close) 
results as the explicit calculations of Ref. [Hj for potential 
widths L > lOl. 
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